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ABSTRACT 

The disruption of a binary star by a massive black hole (MBH) typically leads to the 
capture of one component around the MBH and the ejection of its companion at a 
high velocity, possibly producing a hypervelocity star. The high fraction of observed 
triples (~ 10 % for F/G/K stars and 50 % for OB stars) give rise to the possibility 
of the disruption of triples by a MBH. Here we study this scenario, and use direct 
TV-body integrations to follow the orbits of thousands of triples, during and following 
their disruption by a MBH (of 4x10^ M0, similar to the MBH existing in the Galactic 
Centre; SgrA*). We find that triple disruption can lead to several outcomes and we 
discuss their relative frequency. Beside the ejection/capture of single stars, similar 
to the binary disruption case, the outcomes of triple disruption include the ejection 
of hypervelocity binaries; capture of binaries around the MBH; collisions between 
two or all of the triple components (with low enough velocities that could lead to 
their merger); and the capture of two or even three stars at close orbits around the 
MBH. The orbits of single stars captured in a single disruption event are found to 
be correlated. The eccentricity of the mutual orbits of captured/ejected binaries is 
typically excited to higher values. Stellar evolution of captured/ejected binaries may 
later result in their coalescence / strong interaction and the formation of hypervelocity 
blue stragglers or merger remnants in orbits around SgrA*. Finally, the capture of 
binaries close to the MBH can replenish and increase the binary frequency near the 
MBH, which is otherwise very low. 

Key words: binaries:close-binaries:general-blue stragglers-black hole physics- 
Galaxy:centre-Galaxy:kinematics and dynamics-stellar dynamics 



1 INTRODUCTION 

First theorized by iHillsl l|l988f ). a close interaction by a bi- 
nary star system and massive black hole (MBH) can pro- 
duce a hypervelocity star (HVS) that has sufficient veloc- 
ity to escape the gravitational pull of the galaxy where 
it orig inated, and leave its c ompanion orbiting near the 
MBH (|Gould fc Quillenll2003l'). Since t he discovery of the 



first HVS in 2005 (|Brown etall |2005') , at least 16 HVSs 



Hirsch et al. 


2005 


Brown et al. 


2007 



Brown et al. 2006a|; iBrown et al.ll20061: 



Brown et al., ,20091 ). There are a number 



of mechanisms that can produce HVSs. The Hills' mecha- 
nism has been well studied and appea rs to be th e most likely 
source of HVSs in the Milky Way (|Yu fc Tremaine .2005 : 
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Ginsburg fc Loebl I2OO6I; iPerets et al] l2007l : iPeretsI l2009al : 



Madigan et al.ll2009l . [2OI0 II. Other scenarios include the in- 
teraction of stars with stellar black holes jO'Learv fc Loebl 
I2OO8I ) , and t he inspiral of an intermediate mass black 
hole ( IMBH) jHansen fc MilosavlieTpl2003l : IYu fc Tremaind 



I2OO3I : iLevinlbOOd : ISesana et al.ll2009l ). It was also suggested 
that a binary star system could be ejected as a hypervelocity 
binary through int e ractions with a n IMBH (|Lu et al.ll2007l : 
ISesana et aLll2009l ). [Pereti (|2009bh (hereafter Perets I) ex- 
panded upon the binary hypervelocity scenario and showed 
that massive hypervelocity binaries could be ejected through 
a triple disruption by a MBH. 

Some HVSs have been found at such distances from 
their likely origin at the Galactic Centre (GC) (see 
iBrown et all 120091 '). that despite their high velocities, the 
main sequence (MS) lifetime of these stars is less than their 
flight time. Perets I suggested that this discrepancy can be 
solved if the stars were originally ejected as binaries, and 
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later coalesced to form single rejuvenated stars (blue strag- 
glers) . 

Stellar dynamics (mostly of young B-stars observed 
close to SgrA*) within the central arcsecond of the GC 
provide overwhelming evidence for the exist ence of a cen- 

4 X IO^ Mq (e .g. 'Schodd et al 



tral MBH with mass 



200S: Ghez et al. 



2003l : lReid fc Brunthaleill2004l : lGhez et al.l 2005: G hez et al 
20091 ). ucing HVSs via 



the Hills' mechanism, the tidal breakup of a binary could 
lead to the capture of a star in a close orbit around 
the MBH (|Gould fc QuillenI |2003| ). possibly explaining 
the origin of young B-stars around SgrA* (P erots et al. 
'2007': Madiga n et al.[|2009l:|Peret^ '2009a': Perets fc Fabrvckvl 



2009.: 



Madigan et al.l I2OI0I : IPerets fc Gualandrid I2OI0I ) 



Ginsburg fc Loebl (|2007l ~ showed that the tidal breakup of 



a stellar binary by the MBH can also lead to a collision 
that ends with coalescence. Such mer gers may be a source 
of young stars for th e Galactic Centre (|Antonini et al.|[2010l : 
lAntonini erailbOld ). 

The tidal binary break-up scenario provides a wide va- 
riety of outcomes, which were studied at length. Here we 
expand the study to the disruption of triple stars, which 
given the observed high triple f requency (^ 10 % for F/G/K 
stars and ~ 50 % for OB stars; iRaghavari et al.ll2010l : lEvansl 
I2OIII ). are highly likely to occur. Th is possibility w as dis- 
cussed in general terms analytically (|Lu et al.ll2007l : Paper 
I). We examine the results of over 2000 A'^-body simulations 
of triple disruptions by a MBH, enabling us to study their 
various outcomes quantitatively, and compare them with bi- 
nary disruptions. We discuss these outcomes and their rel- 
ative frequency, as well as their long term evolution. Note 
however, that given the large phase space of possible initial 
conditions, our simulations are not exhaustive, and cover 
only a limited range of conditions. 

In section 2 we describe the codes and simulation pa- 
rameters we used. In section 3 we discuss the triple disrup- 
tions scenario. In section 4 we discuss the formation of a 
blue straggler (BS) and consequently the formation of a hy- 
pervelocity binary (hereafter HVB) from a triple disruption. 



2 COMPUTATIONAL METHOD 

In our s tudy wo have used the N-body code written by 
lAarsethl (,1999 ), which is freely available for download B We 
adopted a small value of 10~* for the accuracy parameter rj, 
which determines the integration step through the relation 
dt = y/ rjF /{d^F/ df^ ) where dt is the timestep and F is the 
force. The softening parameter, eps2, which is used to create 
the softened point-mass potential, was set to zero. We treat 
the stars as point particles and ignore tidal and general rela- 
tivistic effects on their orbits, since these effects are small at 
the distance (~ lOAU; over the short dynamical timescale 
of the encounter) where the triple or binary system is tidally 
disrupted by the MBH. 

We set the mass of the MBH to M = 4x lO^Af©. For our 
initial run, we set the mass of each star to m = AMq within 
the triple system; comparable to the currently observable 
hypervelocity stars. For the next set of simulations, we used 



^ http://www.ast. cam. ac.uk/^sverre/web/pages/nbody.htm 



the same parameters as that of the triple systems, however 
we used a binary system with mass mi — 8Mq and m2 = 
4M0 , equivalent to the combined mass of the triple system. 
Our initial simulations were planar, however we also ran 
over 500 simulations with an inclination angle of 60 degrees. 
Most of our simulations assumed initial orbital eccentricities 
of e = 0. For comparison, we ran a set of over 500 simulations 
with e = 0.4. All runs start with the centre of the circular 
triple system located 2000 AU (= 10~^pc) away from the 
MBH along the positive y-axis. This radius is larger than 
the binary size or the distance of closest approach necessary 
to obtain the relevant ejection velocity of HVSs, making 
the simulated orbits nearly parabolic. We used the same 
initial distance for all runs to make the comparison among 
them easier to interpret as we varied the distance of closest 
approach to the MBH or the relative positions of the stars 
within the system. 

For most of our simulations involving triple systems, 
we chose the inner binary to have semimajor axis of length 
flo — 0.05, 0.1, 0.15, and 0.2 AU, and left the second semima- 
jor axis at 5ao in order to keep the binary as tight as possible. 
However, we did run simulations with larger ao as well as 
simulations with semimajor axis length greater than 5ao to 
see how this increase affects the disrupted systems (see Fig- 
ure [2] and Table |3]). Such a distribution of semimajor axis 
lengths allows us to test the range of hypevelocity stars and 
hypervelocity bina ry producti o n. Wider binaries giv e lower 
ejection velocities |Hillsl Il99ll : iBromlev et al. I I2OO6I ). Much 
tighter binaries would not be easily disrupted by the black 
hole, or may coalesce to form a single star before interacting 
with the MBH. The size of a MS star of a few solar masses 
is ~ O.OIAU, and so binaries tighter than ~ 0.02AU are 
precluded because the two stars will shortly merge. 

In the Galactic disk, ~ 0.75 of all stellar sys- 
tems with O/B and ~ 0.3 of systems with F/G/K pri- 
mar ies, are binaries or higher multiplicity system s (see 
e.g. iRaghavan eralll2010l : iKobulnickv fc Frveij BoOtI ') . with 
roughly equal probability per log arithmic interval of sepa- 
rations, dP/dln{a) = const (e.g. lAbtllTgsl : iHeacoxl Il998l : 
Larson 2 003i ) . F/G/K binaries have a log- normal distribu- 
tion of periods (|Raghavan et al.ll2010l ). For such a distribu- 
tion, ~ 0.3 — 0.4 of all OB binaries have < 1 AU separations 
(see discussion of binary and triple fractions in Paper I), 
where F/G/K binaries are likely to only have ~ 0.05 of the 
systems in this range. 

The initial phase of the stellar orbi t plays a crucial role 
in the outcome l|Ginsburg fc Loebll2006h . Therefore, we sam- 
pled cases with initial phase values of — 180° in increments 
of 15°. As initial conditions, we gave the binary system no 
radial velocity but a tangential velocity with an amplitude 
such that our effective impact parameter is ~ 5 — 10 AU. 
We expect no HVSs to be produced at substantially larger 
impact parameters (Ginsburg fc L oeb 2006). 

We also examine cases for triple systems with tight in- 
ner binaries (0.05 and 0.1 AU) and masses in the range 
1,2,3,5 and 6Mq per star. Note that a triple system of 
Mtripie ^ 3Mq with a close outer binary is far less likely 
than one where the primary is of hig her mass m ~ 5Mq 
(iTokovinin et al ] |2006l : lTokovininll2008l ) . Upon completion of 
all our simulations, we evolved the H VBs using the binary- 
star evolution (BSE) code written bv I Hurley et al.l l|2002h . 
and compared the results with the evolution of a single star 
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of similar mass usi ng the single-star evolution (SSE) code 
jHurlev et al.ll2000l ). Our results show that a HVB can coa- 
lesce well after its ejection, and be observable as a MS star 
even after 4 — 5 times the lifetime of a single star of similar 
total mass (m ~ 7Mq). Such a HVB could therefore be ob- 
served at far larger distances than accessible to a single star 
with the same total mass and velocity, which would evolve 
off the MS. 



3 THE TRIPLE DISRUPTION SCENARIO 

Given a binary system with stars of equal mass m separated 
by a distance a and a MBH of mass Af S> m at a distance b 
from the binary, tidal disruption would occur if fe < bt where 



m 



(1) 



The distance of closest approach in the initial plunge of the 
binary towards the MBH can be obtained by angular mo- 
mentum conservation from its initial transverse speed v± at 
its initial distance from the MBH, d, 



vid - 



1/2 



(2) 



The binary will be tidally disrupted if its initial transverse 
speed is lower than some critical value. 



{GMa 



,1/2 



d 



M 



1/6 



1/2 



— =10" 



1/6 J 



km s 



(3) 

where a_i = (a/0.1 AU), da. 3 = (d/2000 AU), mo.s = 
(m/SMo), and we have adopted M = 4 x 10^ M©. For 
v± < iix.crit, one of the stars receives sufficient kinetic en- 
ergy to become unbound, while the second star is kicked into 
a tighter orbit around the MBH. The ejection speed, v^^ , of 
the unbound star can be obtained by considering the change 
in its kinetic energy ~ vSv as it acquires a velocity shift of 
order the binary orbital speed Sv ~ ^jGrnfa during the 
disruption process of the binary at a distance ~ ht from the 
MBH when the binary centre-of-m ass speed is ii ~ \J GM/bt 
('Hills" 1988': ' Yu fc Tremainell2 003''). At later times, the binary 
stars separate and move independently relative to the MBH, 
each with its own orbital energy. For v < u_L,crit, we therefore 
expect 



1/2 



1/2 



1.7 X 10'' m, 



1/3, 



-1/2 
-1 



km s 



(4) 



For a triple system, the possible outcomes are (i) the 
stars are disrupted and separately orbit the MBH; (ii) a 
HVS is produced and the remaining binary orbits the MBH 
(the binary might then be disrupted in such a way that a 
second HVS is produced); (iii) a HVB is produced and the 
remaining star orbits the MBH; (iv) a collision occurs be- 
tween either two or all three stars; (v) the stars are disrupted 
but a close binary remains around the MBH. During a col- 
lision, if the impact velocity is less than Vesc, the stars will 



Outcome For 3 Stars 


P(a 


^ 0.2 AU) 


HVS -1- 2 Capture 




0.26±0.03 


HVS + Binary Capture 




0.02±0.01 


2 HVS -1- 1 Capture 




0.04±0.01 


HVB 




0.05±0.01 


Collision 




0.04±0.01 


Binary Capture (no HVS) 




0.11±0.02 


3 Capture 




0.48±0.05 


Outcome For 2 Stars 


P(a 


^ 0.2 AU) 


4M0 HVS 




0.20±0.03 


SMq HVS 




0.13±0.02 


Collision 




O.OlitO.Ol 


2 Capture 




0.54±0.05 


Binary Capture 




0.12±0.02 



Table 1. The first section of the table shows the probability of 
various outcomes for a triple system composed of 4Mq stars upon 
disruption by the MBH. Here the binaries are all of semimajor 
axis a ^ 0.2 AU. The third star is at a distance of 5a for all cases. 
Each row shows the values obtained from our initial simulations 
with their corresponding Poisson errors. The first row shows the 
probability that a single HVS is produced, and the remaining 
stars orbit the MBH independently. The second row shows the 
probability that a HVS is produced, and the remaining stars or- 
bit the MBH as a binary. The third row shows the probability 
that 2 HVSs are produced, and consequently the remaining star 
orbits the MBH. The fourth row shows the probability of produc- 
ing a HVB. The fifth row shows the probability of a collision for 
either 2 or 3 stars (note that a collision here does not necessar- 
ily imply coalescence). The sixth row shows the probability for 
a binary to be captured. In this case there is no HVS produced, 
and the remaining star orbits the MBH in an independent orbit. 
The seventh row shows the probability of the MBH capturing all 
3 stars. In this scenario the final orbit of each star is independent 
of the others. The second section of the table shows the proba- 
bility of various outcomes for a binary system composed of 4Mq 
and 8Mq stars. The binary separation is similar to the binary 
separation for the triple scenario, where a ^ 0.2 AU. The first 
row shows the probability of producing a 4Mq HVS, while the 
second row shows the probability of producing a 8Mq HVS. The 
third row shows the probability of a collision. The fourth row 
shows the probability that both stars are captured by the MBH 
in independent orbits. The fifth row shows the probability that 
both stars are captured by the MBH as a binary. 



likely coalesce. If a close binary has semimajor axis length 
such that mass transfer can occur, then a possible outcome 
is a BS. Scenarios ii-iv are illustrated in Figure [l] Table [T] 
summarizes the statistical results from our initial runs using 
stars of 4Mq with a ^ 0.2 AU. 



3.1 Comparison With A Binary System 

In order to verify whether the initial disruption of a triple 
system behaves differently than a binary disruption, we ran 
simulations of binary disruptions by the MBH with the same 
parameters as those in our triple systems. The innermost 
star was given a mass of 8Mq , the companion 4Mq , and all 
other properties (semimajor axis, inclination, impact param- 
eter) remained the same. Although a binary system can not 
produce a HVB, most of the other outcomes are statistically 
similar to those of the triple systems. Our results are sum- 
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Figure 1. The three images show the evolution of a triple system after being disrupted by the MBH. In each scenario the system is 
initially located 2000 AU along the positive y-axis. The separation between the inner binary is a = 0.1 AU, and between the inner binary 
and third star is 0.5 AU. The stars of the inner binary are denoted by the red and blue colors respectively, and the outer binary by the 
color green. Each star has a mass of AMq. In the first case (left image) the triple system produced a HVS upon being disrupted and the 
ejected star was the outer binary (the HVS's path is denoted by the solid green line). The remaining binary orbited the MBH until it 
was also disrupted and a second HVS produced. The path of this second HVS is denoted by the solid red line. The remaining star (blue) 
is in a highly eccentric orbit around the MBH. The middle image shows the production of a hypervelocity binary (HVB). In this case, 
the outer star remains in a highly eccentric orbit around the MBH, while the inner binary (denoted by the magenta star) is ejected. In 
the right image the triple system is perturbed by the MBH, and before completing an orbit all three stars collide with v < Vesc where 
the escape speed for a star of a few solar mass is ~ 500fcms~^. It is likely that in such a scenario the three stars will coalesce. 



marized in Table [T] We conclude that the initial disruption 
of a tight triple system or binary system behave similarly. 

3.2 Inclination and Eccentricity 

Most of our simulations were run with planar configurations 
(relative inclination of the binary and its trajectory around 
the MBH of i = 0°), and with initial binary orbital eccen- 
tricity of e = 0. In order to explore the effects of the inner 
binary properties on the outcomes, we have run additional 
simulations with mutual binary inclination of i = 60° (on 
circular orbits; e = 0) as well as simulations of eccentric bi- 
naries, e — 0.4 (with planar orbits; i = 0°) while keeping 
other parameters unvaried. The orbital inclination of the 
captured star s is determined pri marily from their initial in- 
clination (see IZhang et al.l I2OI0I for a discussion on binary 
HVSs, the case is similar for triple disruptions). Our results 
are summarized in Tabled While roughly similar to our re- 
sults from Table [T] we find that for a high inclination and 
higher eccentricity, more HVSs are produced. In addition, 
higher inclination leads to more collisions. 



4 THE OUTCOMES OF TRIPLE 
DISRUPTIONS 

In the following, we discuss the possible outcomes of triple 
disruptions and their potential observable products. 

4.1 HVB 

The Hills' mechanism which is responsible for the produc- 
tion of HVSs can also produce a HVB through the disruption 
of a triple system. In order to produce a HVB the triple sys- 
tem will need to have Mtripie ~ 6 — 12Mq and atin < 1 AU 



outcome 


P(i = 60°) 


P(e = 0.4) 


HVS -1- 2 Capture 


0.49±0.05 


0.38±0.04 


HVS + Binary Capture 


0±0 


OitO 


2 HVS -1- 1 Capture 


0.08±0.02 


0.02±0.01 


HVB 


0.02±0.01 


O.OSitO.Ol 


Collision 


0.07±0.02 


0.02±0.01 


Binary Capture (no HVS) 


0.04±0.01 


0.10±0.02 


3 Capture 


0.3±0.04 


0.43±0.05 



Table 2. Probability of various outcomes for a triple system com- 
posed of 4Mq with orbital inclination 60 degrees (left column), 
and for a system with initial eccentricity e = 0.4 (right column) 
upon disruption by the MBH. For both scenarios the semimajor 
axis of the third star was 5a where a is the semimajor axis of the 
inner binary with a ^ 0.2AU for all simulations. Statistically, the 
results are comparable to Table [T] with some notable differences. 
In both cases we have more HVSs produced, and less complete 
captures of all stars. Also for the first case we have more collisions 
and less binary capture rates which are due to Kozai-like effects. 
Both cases are meant to illustrate possible situations, and are not 
conclusive in and of themselves. 

(Perets I). Although the fraction of s uch binaries is not well 
known, it is estimated to be ~ 0.05 (|Feke]||l98ll : Perets I). 
A HVB produced this way will be a close or contact binary. 
Consequently, it is possible that most of the observed con- 
tact bina ries were produced through evolution in triple star 
syste ms (jPribulla fc Rucinskil I2OO6I : iFabrvckv &: Tremainel 
l2007f ). Our simulations confirm that all HVBs are close or 
contact binaries with orbital period of ~ 2 days. S uch close 
binar ies are necessary for the formation of BSs (ILu et al.l 
I2OO7D . Furthermore, there is evidence that sug gests nearly 
all ob served BSs are formed in binary systems (|Leigh et al.l 
l2007l ). 

In our simulations we primarily concentrated on tight 
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triple systems (ao ^ 0.2 AU) with stars of mass 4M0, how- 
ever we also ran simulations to see what affect mass and 
semimajor axis length has on the HVB and HVS produc- 
tion rates. We found that an inner semimajor axis length of 
a = 0.2 AU or less is necessary for the production of HVBs, 
so long as the outer semimajor axis is 5a (see Table [3|. Sim- 
ilarly, the rate of HVS production drops dramatically once 
the semimajor axis of the inner binary is > 0.2 AU. Keeping 
the semimajor axes ratio at 5, and the inner binary tight, we 
found that increasing the mass of the system also increases 

the likelihood of a HVB. 

For some bound HVSs fseelBrown et al.ir2007l ). and also 
HVS HE 0437-5439 (|Edelmann et alj 120051 ) there exists a 
discrepancy between their kinematics and apparent ages. A 
number of theories have been proposed to explain the fact 
that their lifetime is longer than expected. Our simulations 
show that although rare (see Table [1]), HVBs can be pro- 
duced via the disruption of triple systems with tight bina- 
ries. Taking the parameters (orbital period and eccentricity) 
from our A^-body simulations and evolving the HVB using 
BSE, we found that the two 4M0 stars will form a common 
envelope, and via mass transfer one star will reach a mass 
of ~ 7Mq. The second star typically goes supernova (mass- 
less), and the remaining star on average stays on the MS for 
~ 200 Myr. Typically, a star of 7Mq will live for no longer 
than ~ 50 Myr, therefore the HVB evolves into a runaway 
BS. Our BSE simulations show that short periods (^ 10 
days) but high eccentricity (e ^ 0.6) will produce a star of 
~ 7Mq with lifetime ^ 50 Myr. The average orbital period 
of our HVBs is ~ 2 days, and the eccentricity is typically 
low, with an average value of ~ 0.3. Such scenarios allow 
for rejuvenation and can explain the unexpected longevity 
of such HVSs. 



4.2 Close Binary Around SgrA* 

The disruption of a triple system is three times more likely 
to produce a captured binary than a HVB (see Table [T|. 
We define a captured binary as a binary with semimajor 
axis ^ 2 AU orbiting the MBH. Our simulations show that 
if two stars remain bound together after their initial tidal 
interaction with the MBH, then they almost always have 
semimajor axis ^ 2. Otherwise the system is disrupted and 
the stars orbit the MBH independently and mass transfer 
is not possible. The orbital parameters of such a binary are 
similar to those of a HVB, and consequently the evolution 
will in most cases be identical to that of a HVB. Our average 
period for a close binary was ~ 10 days, and the eccentric- 
ity ~ 0.4. Such values lead to an average lifespan of ~ 190 
Myr for a BS of mass ~ 7Mq. As in the case for a HVB, 
the eccentricity greatly affects the lifetime of the BS. Tidal 
stripping may also affect the evolution of the binary, and 
it has been suggested that the S-stars orMting SgrA* are 
tidally stripped AGB stars jPavies fc Kind[2005h . however 
there is no confirmation of any stripped stars at the Galactic 
Centre, and current evidence indi cates that S2 is a ge nuine 
early B star of mass 14 - 2OM0 (|Martins et al.l 120081 ). Al- 
though not all of the stars orbiting SgrA* could come from 
triple disruptions, it is possible that some of the stars are 
blue stragglers that evolved from triple disruptions. 



4.3 Collisions 

Direct collisions are one of the primary mechanism for BS 
formation, and in den se environments colli sions often involve 
more than two stars (|Fregeau et al.ll2004l ). Under some cir- 
cumstances a triple system is disrupted in such a way that 
either two or all three stars collide. Assuming that a binary 
is disrupted at a random angle and ignoring gravitational 
focusing, the probability for a collision is 



2'KRse 



(5) 



where i?* is the radius of the star (assuming radii are ap- 
proximately equal) and Rsep is the average separation of 
the stars. The stars will merge if the relative velocity of im- 
pact is less than the escape velocity from the surface of the 
star 500kms~^). In our simulations around 15 per cent 
of all collisions have impact velocities low enough to allow 
two stars to coalesce, and roughly 4 per cent allow all three 
stars to coalesce. A simple estimation for the impact velocity 
comes from conservation of energy. 



E : 



1 17111712 



Gmi7n2 



constant 



2 mi + m2 r 
which yields the relative impact velocity for two stars 



Vf = [2G(mi +m2){- 



1 



1/2 



(6) 



(7) 



The equations are nearly the same for three stars, with 
an added ma. Equation (7) does not take tidal forces into 
account nor gravitational focusing. The rate of HVS pro- 
ductio n in the Milky Way is ~ 10~^yr~^ (jBrown et al.l 
l2006bl). and the total coUisional rate is ~ 10 per cent 
( Ginsburg fc Loebll2007l ). Given that mergers occur in ~ 20 
per cent of all collisions, the expected production rate of 
such S-stars is once very 5 Myr. The lifetime of a massive 
star of mass ~ 14M0 is roughly 10 Myr. Therefore merg- 
ers can not account for a majority of the S-stars, but may 
account for some and can help explain the longer than ex- 
pected lifetimes of these stars. 



5 EVOLUTION OF CAPTURED AND 
EJECTED BINARIES 

Following the triple disruption, captured binaries could 
evolve through short (orbital) and long (secular) term dy- 
namical evolution induced by the MBH on other stars in 
their environment. They may also evolve later through bi- 
nary stellar evolution. HVBs may similarly evolve through 
binary stellar evolution, but their mutual orbits are unlikely 
to be affected by interactions with other stars which are 
scarce once a HVB is ejected from the GC. 



5.1 Dynamical evolution of captured binaries 

The dynamical channels could lead to the disruption of the 
binaries or to their merger /collision. In many cases a triple 
disruption is followed by the capture of a binary, which it- 
self is disrupted shortly after to leave two bound components 
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around the MBH which are unbound to each other (see Ta- 
ble [1] and Table [2]). This possibility leads to both binary 
components orbiting the MBH in individual, but correlated 
orbits. Over longer timescales, relaxation processes such as 
resonant relaxation could change the orbits of the captured 
stars such that t heir inclination an d eccentricity no longer 
seem correlated jPerets et aLlbOO^ ). 

In other cases (see 31}, following the triple disruption, 
a binary is captured but collides to form a merger prod- 
uct. When binaries are captured with high inclination be- 
tween their mutual orbits and the orbit around the MBH, 
they are more likely to collide, as secular Kozai-like effects 
can excite their eccentricity inducing the merger (Paper I; 
lAntonini et al.ll201(]| : Untonini et al.ll2010l '). Indeed, we find 
that triple disruptions of inclined triples (i.e. leading to the 
capture of inclined binaries) show a higher rate of collisions 
at the expense of the binary capture rate, which is decreased 
(see Table [2]). The collision/merger products could poten- 
tially be observed as atypical, more massive S-stars, possibly 
with high rotation velocities. 
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0.2 0.4 0.6 0.8 1 



0.2 0.4 0.6 0.8 1 



Close Binary 



0.2 0.4 0.6 0.8 1 



Binary Separation (AU) 



5.2 Stellar evolution of captured and 
hypervelocity binaries 

Captured binaries which did not collide shortly after their 
capture, as well as ejected HVBs, may later evolve through 
stellar evolution. Such evolution may lead to their merger 
through mass transfer and common envelope evolution, or to 
the production of compact binaries in the GC; or in the case 
of HVBs, to hypervelocity merger products and compact 

binaries. 

First discovered bv lSandaed l| 19531 ). a BS is a MS star 
with an apparently abnormally long lifetime. BSs are mas- 
sive stars that lie above the turnoff region in the blue end 
of the color-magnitude diagram. If BSs were normal MS 
stars, they would have evolved away from the MS. There 
are a number of theories proposed to explain the unexpected 
longevity of BSs. One theory involves stellar collisions and 
coalescence, another mass transfer from a binary compan- 
ion, and a third is from induc ed merger /co llision through 



evolution in triple systems (see IStrvkeill 19931 for a review of 



the two former theories and IPerets et al.ll2009l : Paper I; and 
lAntonini fc Peret3l201ll . regarding the evolution in triples). 

Using the BSE binary stellar evolution code by Hurley 
ijHurlev et al.ll2002h . we evolved two stars of mi = = 
4Mq with period P ~ 10 days. Such a system evolves into a 
7Mq star with a companion of m < IMq . The lifetime of the 
evolved 7Mq star is ~ 150-220 Myr. The lifetime of a single 
a single star with mass 7Mq is ~ 40 Myr. On average, the 
lifetime of the evolved star will be ~ 4 — 5 times longer than 
a single star of comparable mass. Hypervelocity binaries, as 
well as similar products captures around the MBH could 
thus form BSs. However, the latter binaries do not evolve 
in isolation, as they may be affected b y encounters with 
other stars near the MBH (|Peretj|2009al ). or be affected by 
the MBH itself (Paper I). Further study of binary stellar 
evolution in the GC environment is required to predict the 
fate of such binaries. 



Figure 2. This cumulative plot shows the fraction of HVSs, 
HVBs, collisions, and tight binaries produced from the disruption 
of a triple system by the MBH given an initial binary separation 
(AU). These results are based upon our initial simulations with 
each star of mass 4Mq. The production of HVSs, HVBs, and close 
binaries all drop exponentially with increasing binary separation. 
In particular, the production of HVBs occurs only when the ini- 
tial binary separation is very tight. The rate of collisions stays 
fairly constant with increasing binary separation. Note that for 
these simulations the initial separation of the third star is 5 times 
the binary separation. 



6 CONCLUSIONS 

In our study of triple disruptions by a MBH, we mapped the 
various possible outcomes, including the capture and ejec- 
tion of single and binary stars (including HVSs and HVBs), 
collisions between two or all of the stars, and the pertur- 
bations of the orbital parameters of captured/ejected bi- 
naries. We also followed the long term stellar evolution of 
captured/ejected binaries which could lead to their merger. 

Our N-body simulations indicate that if a triple sys- 
tem of inner binary separation 0.05-0.2 AU with the third 
star at five times that distance approaches SgrA* to within 
< 10 AU, the subsequent tidal disruption can lead to the 
production of a HVB. Such HVBs are always tight, and 
stellar evolution code shows that they evolve into massive 
BSs {M final ~ 7Mq). It is possible that some of the bound 
HVSs are rejuvenated stars. However, it has not yet been 
confirmed that all these stars are true B-type stars and 
not horizontal branch stars. HVS HE 0437-5439, has been 
identified as a true young B2 III-IV halo star with mass 
~ 9Mq (|Bonanos et al.1 l2008l : IPrzvbilla et~a l. 2008). Sinc e 
the star likely originated from the GC ( Brown et alll2O10l ). 
it is quite possible that this HVS is a rejuvenated star, orig- 
inally ejected as a binary from the GC (see Perets I). 

In our simulations involving triples, the inner binary 
always has a mass twice that of the third companion. As 
expected, when a binary is ejected it typically has lower 



ejection velocity than compared with that of the companion 
(|Bromlev et al.l [20061 ). The ejection of a HVB in our simu- 
lations is therefore less likely to occur than a HVS ejection. 
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outcome 


ao/fli = 5 


a-o/ai = 7.5 


ao/a.i = 8.75 


ao/ai = 10 


HVS 


0.62±0.1 


0.56±0.1 


0.35±0.08 


0.36±0.08 


HVB 


0.16±0.0.05 


0.02±0.02 


0.04±0.03 


0±0 


Collision 


0.02±0.02 


0.02±0.02 


0.05±0.03 


0.02±0.02 


Close Binary 


0.2±0.06 


0.27±0.07 


0.27±0.07 


D.33±0.08 



Table 3. Probability of a HVS (first row), HVB (second row), collision (third row), and formation of a close binary (fourth row), from 
the disruption of a triple system of AMq stars with varying outer semimajor axis length. The semimajor axis of the inner binary is 
fixed at O.ObAU . This particular inner binary is very tight and thus has a relatively high probability of producing a HVB. The first 
column shows the probabilities of an outcome with an outer to inner semimajor axis ratio of 5. Note that this is the smallest stable 
outer semimajor axis length. The ratio of the semimajor axes in the second column is 1.5 times that of the first column. The ratios of 
the semimajor axes in the third and fourth columns are 1.75 and 2 times that of the first column respectively. All associated errors are 
Poisson. We can see that the production of HVSs drops as we get into larger ratios. However, at a certain distance the production of 
HVSs becomes solely due to the tight inner binary. The production of HVBs drops to ~ for larger outer semimajor axes. The collision 
rate remains constant within error, as does the production of close binary stars. We ran simulations with similar semimajor axis ratios 
but larger inner semimajor axis and results were similar to within the expected error. 



even though the inner binary and the third companion have 
a comparable chance of being ejected following the triple dis- 
ruption. Nevertheless, the HVB ejection rate could extend 
to a non- negligible fraction of the HVSs (up to 1/4 of all 
HVSs in some of the cases we studied) . A higher mass ratio 
between the third companion and the inner binary (as well 
Eis overall more massive triples, even with the same mass 
ratio) would lead to a higher HVB ejection rate. The overall 
behavior is consistent with the assumptions used in previ- 
ous studies of triple disruption (Perets I) , and confirms their 
validity. One observed difference, however, is the possibility 
of a collision between the inner binary components, not ac- 
counted for in the previous studies of triple disruptions. This 
interesting outcome reduces the HVB ejection rate, but only 
by a small fraction. 

As mentioned above, in most cases a HVB is not formed 
upon the disruption of the triple system. However, even if a 
HVB is not formed, there is a fair chance that the resulting 
disruption will produce a tight binary around SgrA* (see 
Table [1} • Barring tidal stripping or close encounters with 
massive systems, a tight binary orbiting SgrA* will evolve 
into a BS in the same manner as a HVB. It is possible that 
some of the massive S-stars are rejuvenated stars. It is also 
possible for a triple system to be disrupted so that two or all 
three stars collide. In some instances these collisions occur 
with V < Vesc and thus the stars can coalesce and form a 
BS. If all three stars coalesce, this mechanism may explain 
some of the most massive stars we see around SgrA*. 

In each scenario discussed in section [5] any star 
that remains orbiting SgrA* is always on a highly ec- 
centric orbit with eccentricities ranging from e = 0.9 
to almost unity. Such high eccentricities are inconsistent 
with the observe d orbital parameters of the S-stars (e.g. 
iGhez et al.ll2009l ). However, resonant relaxation processes 
could change the eccentricity distribution over time, mak- 
ing the capture s cenario a viable model for the origin of 



the GC S-stars llPerets et al.l [20071: iMadigan et al.l 
Peretsll2009al: iPerets &i FabrvckvlbOOg -] iMadigan et al.l l 
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Perets fc GualandrisI |2010| ). We note that the initial 
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ec- 
centricity distribution of captured stars, following the bi- 
nary/triple disruption in our simulations, suggests a lower 
eccentricity (e = 0.94) than previously assumed for the ini- 
tial conditio ns in simulations of the l ong term evolution o f 
the S-stars (jPerets fc Fabrvckvl[2009l : IMadigan et"al]|2010^ . 



Lower initial eccentricities are likely to provide better consis- 
tency between the observations and the simulation results. 

Finally, this work covered only a small part of the possi- 
ble phase space for triple disruption, but allowed us to map 
out their potential outcomes, and validate previous assump- 
tions regarding triple disruption using detailed few-body 
simulations. Future extensions of this work may examine the 
disruption of triple systems of varying masses, inclinations, 
and eccentricities. 
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